Morphological changes of a superconducting phase in a mixed state with a normal 

current: a multiple scale analysis 
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The mixed state of superconducting (SC) and normal (N) phases in one dimensional systems are 
characterized by several phase slips and localization of the order parameter of the SC phase. The 
phenomenon is explained on the basis of a complex Ginzburg Landau (CGL) model. We present a 
simple analysis of the system on multiple scales to capture localization and phase slips when phases 
coexist. 
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Below the critical temperature a one dimensional su- 
perconducting channel is driven from superconducting to 
normal phase by increasing the current passing through 
it. The phase transition is current driven and the role 
of excitations to drive such a transition has been ex- 
plained on the basis of a (CGL) model by Langer and 
Ambegaokar which was later improved by McCumber 
and Halperin (LAMH) [1, . Although LAMH approach 
was phenomenological, similar models have been derived 
on a microscopic basis 0]. The LAMH theory is reck- 
oned to have adequately captured the essential thermo- 
dynamic features of this current driven transition. How- 
ever, this classic problem is often revisited from different 
view points. A recent work that characterizes the pe- 
riodic appearance of phase slip center (PSC) in the su- 
perconducting phase by collision of real eigenvalues at a 
Hopf bifurcation as a consequence of the PT (Parity and 
Time reversal) symmetry of the system is quite inter- 
esting A closed expression for the temperature and 
current dependent free energy barrier of clean supercon- 
ducting micro channels has also been derived recently Q . 

In the present paper we will investigate the effect of 
a small normal (Ohmic) current on a superconducting 
phase in a region where the phases can co-exist. We will 
employ multiple scale perturbation to find the form of 
slowly varying amplitude of the SC order parameter. Al- 
though the analysis is extremely simple yet reveals an ex- 
act form of localization of the SC phase and captures the 
periodic phase slips. We can also draw a few important 
conclusions from this very simple analysis. (1) The non- 
linearity of the CGL does not play any role whatsoever on 
the generation of the localization and the dynamic phase 
slip centers (PSC). (2) Phase slips are a consequence of 
the spatial variation of the amplitude of the order pa- 
rameter than a direct consequence of the existence of the 
Ohmic current. (3) The localization of the SC order is 
a direct consequence of the existence of the Ohmic cur- 
rent. (4) There are indications that the local nature of 
the phase slips are a consequence of PT symmetry of the 
system in the absence of which it either have not occurred 
or have been smeared over space. 



The CGL model for ID Superconductor reads as 

u (-0* -I- zmV') = i^xx + (1 - 



(1) 



where if) is the superconducting order parameter. The 
system being ID, x is the distance along the wire, /i is 
the electrochemical potential which characterizes the N 
phase of the system. The j is the current and parameter 
u is the ratio of order parameter relaxation time to the 
current relaxation time. The typical value of u is 12 in the 
strong depairing limit whereas 5.79 for weak depairing 
^6]. We will treat u as a constant here. The suffix t and 
X oi indicates corresponding derivatives of it. 

Eq.l has two stationary solutions. (1) ip = 0, = —xj 
(for constant j) which is the normal state and (2) ~ 
Aexp{iqx), — 1 ~ , j = A^q which is the supercon- 
ducting state when /i = 0. There is a critical current 
jc above which the SC state gives up to the normal state 
and just below jc over some range on the current scale the 
SC phase is metastable and often seen to coexist with the 
normal state. An extensive numerical analysis bvKramer 
et al in this region showed localized SC phase [7|] whose 
profile very much resembles a Gaussian one. The corre- 
sponding profile of the /i is like a straight line passing 
through the center where the SC profile peaks. In what 
follows, we will analyze the system in such a region below 
the jc (say) where two such phases can coexist. 

Since the dominant phase below jc is the SC phase, 
let us consider its order parameter ipse to have the form 
Ipse — A{X, T)[tp° + eip^] where the amplitude of the order 
parameter A{X, r) is a function of slow space and time 
scales X and r respectively. We also take into account an 
additive small variation (0(e)) on the order parameter in 
the form of -0^ which supposedly appears in the presence 
of an equally small (0(e)) presence of a normal current 
(ejAf). We are interested in getting a form of A{X,t) in 
such a situation. Since we perturb the steady SC state, 
the time scale under effective consideration here is the 
slow one only t = er, but we should take multiple length 
scales + eX where the amplitude is the only 

slowly varying function. 
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In the present analysis we will neglect the role of non- 
linearity to enable us apply multiple scale perturbation 
technique as it is normally done near a linear instability 
threshold to derive slow amplitude equations. At this 
stage, let us modify the model Eq.l in such a way that 
does not affect the basic steady state results of the sys- 
tem. We will from now on consider in the place of 
It/'P in the model where A is exactly the constant ampli- 
tude of the steady SC order parameter. This will render 
the model linear without changing the steady state solu- 
tions. So, form now on in what follows we will treat A 
as a constant if not explicitly mentioned to be function 
of slow variables. 

Now, at 0(1) we get 

VL + (1-^^)V° = (2) 

whose solution is the steady superconducting state = 
A{X,T)exp{iqx) where at smaller scales A{X,t) is effec- 
tively constant A. Thus, the above made modification 
of CGL can also be viewed as considering dynamics on 
widely varying scales. At 0(e) we get 



V'l, -t- (1 - ^2)^1 ^ 



dA{X,T) dA{X,T) 
\u i2q- 



dr ^ dX 

- ijNxA{X,T)]exp{iqx), (3) 

where on the left hand side we have absorbed the factor A 
of tp^ within ip'^ since A is effectively treated as constate 
on small scales. The term on the right hand side of the 
above equation is a secular term. So, it has to vanish 
for ip-^ to exist since the linear operator (considering A^ 
in the place of IV'^P) on the left hand side has a zero 
eigenvalue solution. Thus, we get the linear amplitude 
equation 

dA i2q dA iejN 
'dt ~ 



-xA = 



(4) 



u dx u 

The e in the last term appears as a consequence of scale 
change to the normal ones. Consider the transforma- 



tion A = Ae 



So, 
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-2mxAe 



-mx dA 



Putting this in the above amplitude equation we get a 
selection for m as m = which will cancel the terms 
with a factor x in the amplitude equation. The inverse 
characteristic length m is positive for positive ^ and 
causes a Gaussian localization of the SC order parame- 
ter. Its important to note that for constant normal cur- 
rent jAT the electrochemical potential is a straight line 
passing through the origin in space and our analytic ex- 
pression very closely resembles what has been numeri- 
cally obtained by Kramer et al in [3] . The localization of 
the SC phase results as a consequence of the existence of 
nonzero normal current. An interesting point to note is 
that the smaller the normal current the wider is the SC 
phase. The SC state also spreads for larger q. 

Let us concentrate on the rest of the amplitude equa- 
tion. It reads 

dA{x,t) i2qdA{x,t) 



Consider the form of A{x,t) as A{x,t) = a{t)exp{±Xx) 
with A real. So, we readily get A{t) = exp{i^^t) where 
the amplitude oscillates with a frequency uj — 
Since, the amplitude equation is linear, a linear superpo- 
sition of the solutions are always possible and a general 
form of the amplitude can be written as 



A{x,t) 



r Xx iujt 

e e 



e e \e 



(6) 



dt 



dx 



(5) 



This amplitude has a real temporal oscillation at a; = 
in the form of Cos{ujt) due to superposition. Since, the 
amplitude of the two temporally oscillating terms vary 
every where except at the origin in space we will only 
get localized real oscillatory solutions at this point. This 
localized oscillation will cause periodic phase slip at ori- 
gin when the amplitude vanishes. Interestingly enough, 
the period of oscillation is not dependent on the and 
thus phase slips are always there even if we do not con- 
sider the presence of the normal current from the very 
beginning. It only depends on the exponential variation 
of the order parameter amplitude in space. However, 
in the presence of the normal current the localization of 
the amplitude will suppress the exponential divergence 
of the other part at large distances and thus would make 
the solution an acceptable one. Nevertheless, one can 
easily guess that exponential amplitude profiles formed 
locally due to inhomogeneities or other causes can also 
cause periodic phase slips. 

Another important point to note is the joint variation 
of sign in the space and time part of the linear solution 
is actually enabling us to have the localized oscillation. 
If we had uj of both the signs corresponding to the same 
sign of A the superposition would have resulted in a global 
oscillation. This fact indicates that the symmetry of the 
system under joint inversion of space and time is impor- 
tant for generation of localized phase slips which has been 
investigated numerically at length in [j| . The other very 
important point revealed by our analysis is that there 
is no need of the nonlinearity for the generation of the 
localization and phase slip centers. The nonlinearity def- 
initely plays its role on the stability of these structures 
as has been revealed by previous works [H, [3, [1] . 

We will conclude by saying that a very simple multiple 
scale analysis of the of the CGL model for ID supercon- 
ductors analytically captures most important character- 
istics of the system near a region close to the transition 
from SC to normal phase. (1) It reveals the nonlinear- 
ity is not at all responsible for intrinsic generation of 
localization and PSCs of the system. (2) The localiza- 
tion of the SC phase is a direct consequence of the pres- 
ence of the normal current whereas the PSCs are not. 
(3) PSCs will appear even at a zero normal current but 
with suitable local variation of SC amplitude in space (if 
the normal current does not indirectly results from the 
amplitude variation) caused by inhomogeneities or other 
causes. (4) Symmetry under joint space time inversion 



3 



seems to play an important role in the localization of the 
periodic phase slips. 
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